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The conducting and mechanical properties of a metallic nanowire formed at the junction between 
two macroscopic metallic electrodes are investigated. Both two- and three-dimensional wires with a 
W(ide)-N(arrow)-W(ide) geometry are modelled in the free-electron approximation with hard-wall 
boundary conditions. Tunneling and quantum-size effects are treated exactly using the scattering 
matrix formalism. Oscillations of order Ef/Xf in the tensile force are found when the wire is 
stretched to the breaking point, which are synchronized with quantized jumps in the conductance. 
' The force and conductance are shown to be essentially independent of the width of the wide sections 

(electrodes) . The exact results are compared with an adiabatic approximation; the later is found to 
overestimate the effects of tunneling, but still gives qualitatively reasonable results for nanowires of 
length L ^> \f, even for this abrupt geometry. In addition to the force and conductance, the net 
charge of the nanowire is calculated and the effects of screening are included within linear response 
theory. Mesoscopic charge fluctuations of order e are predicted which are strongly correlated with 
' the mesoscopic force fluctuations. The local density of states at the Fermi energy exhibits nontrivial 

behavior which is correlated with fine structure in the force and conductance, showing the importance 
of treating the whole wire as a mesoscopic system rather than treating only the narrow part. 
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PACS numbers: 73.23.Ps, 62.20.Fe, 73.40. Jn 



I. INTRODUCTION 



CO 

G\ ■ Metallic nanowires may be formed at the junction between two metallic electrodes which are pressed together and/or 
pulled apart in a controlled fashion!! a In a wire of nanoscopic cross section, the transverse motion is quantized, leading 
to a finite number of electronic modes below the Fermi energy Ep which can be transmitted through the wire. A 
striking consequence of these discrete modes is the quantization of the wire's conductance at integer multiples of Go = 
2e 2 /h, a phenomenon first observed in two-dimensional (2D) .sejHiiconductor heterostructures,LTEl and subsequently 
studied in detail in three-dimensional (3D) metallic nanowires. ErtM A subtlety inherent in conductance quantization 
experiments is that even for a nearly ideal nanowire, the presence of disorder in the electrodes far from the region of 
interest leads to a suppression of the conductance plateaus below integer values. This suppression is normally taken 
into account by subtracting a phenomenological series resistances*! which allows one to shift the experimentally 
observed plateaus back to integer values. Theoretical histogramai3~t3 exhibit a similar shift towards lower values of 
the conductance, though the pre cise form of the suppression is not equivalent to a simple series resistance due to 
quantum interference effect sJ13'0 These considerations underline the importance of treating the nanowire and the 
electrodes connected to it as a single mesoscopic system, rather than considering the nanowire in isolation. 

The cohesive properties of good metals are also determined to a large extent by the conduction electrons. Thip^ 
one may expect the mode quantization in a nanowire to have a strong effect on its mechanical properties as well.E£l 
In a pioneering experiment published in 1996, Rubio, Agrai't, and Vieira measured simultaneously the force and 
conductance during the formation and rupture of an atomic-scale Au nanowire.cl They observed oscillations in the 
tensile force of order InN under deformation, which were synchronized with jumps of order 2e 2 jh in the, conductance. 
Similar experimental results were obtained independently by Stalder and Diirig.cl In a previous paper JIj it was shown 
that this intriguing behavior can be understood quantitatively using a simple free-electi'on jellium model for a metallic 
nanowire. The theoretical approach introduced by Stafford, Baeriswyl, and BiirkiEj uses the electronic scattering 
matrix to describe the coupling of the nanowire of interest to the macroscopic probes (e.g., STM/AFM) used to 
manipulate it; the correct treatment of this coupling is crucial for calculating the mesoscopic corrections to the bulk 
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electrical and mechanical properties. In the previous paper, the scattering matrix was evaluated using the adiabatic 
and WKB approximations, appropriate for a smooth geometry in which the diameter D of the nanowire varies slowly 
along its symmetry axis z, i.e., (dD/dz) 2 <C 1. The qualitative picture which emerged from the analysis of Ref. 
|l5| is that each quantized mode contributing 2e 2 jh to the conductance of the nanowire also contributes an amount 
of order EpjXp to its cohesive force, where is the De Broglie wavelength of electrons at the Fermi energy Ep. 
(For monovalent metals, Ep/Xp is of order InN.) Under elongation, the cross section of the nanowire narrows, and 
each time a transverse mode is cut off, both the conductance and the cohesive force decrease abruptly. Using an 
elegant argument based on the technique of Ref. |l5|, floppier and Zwerger showed that the leading-order mesoscopic 
correction to the tensile force depends only on the topology of the cross section of the nanowire.113 The scattering 
matrix formalism has also been used to study the effects of impurity scattering on the mechanical properties of 
nanowires .EH I— . 

Independently, van Ruitenbeek, Devoret, Esteve, and Urbinalia considered wires with a special W(ide)-N(arrow)- 
W(ide) geometry, with two wide outer sections, representing the electrodes, and a narrow inner section of constant 
diameter, representing the nanowire (see Fig. |l|). They considered the limit (cf. also Rcfs. |f9|,|o|) where the narrow 
section is sufficiently long that the boundary effects at the junctions of the wide and narrow sections give a negligible 
contribution to the energetics of the nanowire. They also pointed out that screening should be included, and imposed 
a charge neutrality constraint as a first approximation. However, they noted that the assumption of local charge 
neutrality of the nanowire breaks down for very short wires, such as those investigated in the experiments of Refs. 
which are only of order lnm in length. For such short wires, boundary effects are important. 

In the present article, we investigate both 2D and 3D wires with the WNW geometry, treating the boundary effects 
arising from the connection of the nanowire to the electrodes exactly via the scattering matrix approach. Our results 
for the force in the WNW geometry are qualitatively similar to the results for smooth geometries presented in Ref. [D| 
although the presence of sharp corners is shown to lead to a non-generic topological correction to the force for small 
deformations from an ideal wire, which is well understood. The exact results for the WNW geometry are compared 
with an adiabatic approximation, which is found to overestimate the effects of tunneling, but still gives qualitatively 
reasonable results for the conductance and force of nanowires of length L > Af, even for this abrupt geometry. The 
local density of states at the Fermi energy is also calculated, and is shown to exhibit strong modulation on a-Jength 
scale of order A^. This electronic structure of the scattering states gives vivid pictorial support of the notionEj that 
conductance channels act as metallic bonds. 

In addition to the force and conductance, the net charge of the nanowire is calculated, and the effects of screening 
are adressed by a linear response approach. We predict mesoscopic charge fluctuations of order the fundamental charge 
quantum e, which are synchronized with the quantized steps in the wire's conductance, and strongly correlated with 
the mesoscopic force fluctuations. Similar charge fluctuations are predicted in 2D and 3D nanowires; they should 
thus be present in quasi-two-dimensional quantum point contacts exhibiting conductance quantizationTa as well. 
The smallness of the predicted mesoscopic charge imbalance leads us to neglect electron-electron interactions in 
our treatment of the conducting and mechanical properties of the nanowire. Indeed, we find that the mesoscopic 
corrections to the force for wires of length L ^ Xp are large compared to the corrections expected due to screening 
effects, justifying our approach. 

The paper is organized as follows: In Sec. H, we review the scattering matrix formulation of electrical conduction 
and statistical mechanics. General expressions for the force, charge and conductance of a mesoscopic conductor in 



terms of the electronic scattering matrix are derived. In Sec. Ill, the scattering matrix for 2D and 3D wires with 
the WNW geometry is calculated. For simplicity, 3D wires with a square cross section are considered, but it is 
straightforward to use our method for wires of arbitrary cross section if the eigenfunctions and eigenvalues of the 2D 



Schrodinger equation with Dirichlet boundary conditions for this shape are known. In Sec. IV, the results for the 
force, charge, and conductance of a 2D nanowire are presented. The semiclassical approximation to the force and 
charge, and the topological contribution due to sharp corners are discussed. The exact results are compared to an 
adiabatic approximation, and the local density of states at the Fermi energy is calculated. In Sec. [v| the results for 
the force, charge, and conductance of 3D nanowires are presented. The effects of screening are evaluated within linear 
response. The relevance of our results to the experiments of Refs. ^ and [| is discussed. Some general conclusions are 



presented in Sec. VI, and a comparison of the jellium model used here to atomistic descriptions of nanowires based 
on classical molecular dynamics simulationsEiTES is given. 

II. S-MATRIX FORMALISM 

A nanowire connecting two macroscopic electrodes (depicted schematically in Fig. I) is an open quantum mechanical 
system. The Schrodinger equation for such an open system is most naturally formulated as a scattering problem. 
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The fundamental quantity describing the properties of the system is the scattering matrix S(E) connecting incoming 
and outgoing asymptotic states of conduction electrons in the electrodes. (The degrees of freedom corresponding to 
the metallic ions and core electrons will not be treated explicitly, but will be assumed to give rise to a confinement 
potential for the conduction electrons, i.e., to specify the geometry of the system. This should be a reasonable starting 
point to describe simple monovalent metals.) .-The formulation of electrical transport in terms of the scattering matrix 
was developed by LandauercJ and Biittiker^J while the formulation of the statistical mechanics of open quantum 
systems in terms of the scattering matrix was first given by Dashen, Ma, and Bernstein,c3 and was recently revived 
in the context of the persistent current problem by Akkermans et a/.E3 A unified treatment of the electrical and 
mechanical properties_qf metallic nanostructures in terms of the electronic scattering matrix was given by Stafford, 
Baeriswyl, and Biirki.lij In the remainder of this section, we recapitulate the general formalism of Ref. [l^, which will 
serve as the starting point for the present investigation. 

The essential ingredient in the scattering matrix description of mesoscopic systems is that electrons are injected 
into the system from macroscopic reservoirs in internal thermal equilibrium; any perturbation of the reservoirs due 
to the mesoscopic current Sowing from one to another is assumed to be negligible. The energy distribution of the 
electrons injected from reservoir a is thus given by the Fermi distribution function f a {E) = {e~scp[f3(E — + 1} , 
where fi a and /3 = 1/fceT are the electrochemical potential and inverse temperature, respectively, of reservoir a. 

The asymptotic scattering states of conduction electrons for the geometry depicted in Fig. 1 are described by a 
transverse quantum number v and a wave number k which is a function of energy E and v. 



k v {E) = ^{E-e v ), (1) 

where t v is the energy of the transverse modes. If the amplitudes of incoming currents at energy E are given by a 
vector with components a+ (currents from the left side) and a~ (currents from the right side), the outgoing current 
amplitudes b~^, and b~, are given by 

0= s(£ >("-)' (2) 

The scattering matrix for a two-terminal conductor may be decomposed into four submatrices 

S(E) = ( ) , (3) 

where the submatrices Sij(E) contain transmission {i ^ j) and reflection (i = j) amplitudes. Each submatrix Sij(E) 
is a matrix in the incoming and outgoing scattering channels v and v' . 

The electrical conductance is given in terms of the scattering matrix by the well-known formula 

G=2 ~vJ dE ^P Tr { S l2(E)SME)} , (4) 

where the factor of 2 accounts for spin degeneracy. Eq. (^) indicates that G is proportional to the sum over the 
transmission probabilities of electrons incident in a window of width fc^T about the Eermi energy fi. Eq. may be 
modified by electron-electron interactions at finite temperatures, but has been shownEj to hold quite generally in the 
limit T 0. 

The quantity needed to investigate the statistical mechanics of the nanowire i&nrthe electronic density of states, 
D(E). The density of states can be expressed in terms of the scattering matrix ascSEa 

D(E) = _L^ t(E) ^_ h .,). ,„ 

This expression holds for an arbitrary interacting gas of particlesS Given the density of states, the grand partition 
function may be evaluated, and the thermal expectation values of all observables may be calculated. 
The expectation value of the total electronic charge on the nanowire is given by 

(Q-) = -eJdEf(E)D(E). (6) 

Integrating by parts, and taking the limit T — > 0, one obtains the simplified expression 
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(Q-) = --Im{]ndetS(E F )}. (7) 

7T 

The overall phase of the scattering matrix depends on the exact choice of the asymptotic states (their phase) . Therefore 
the phase relation chosen between the amplitudes and bf, is in principle free. Different choices correspond to the 
inclusion of various amounts of the constant asymptotic charge density in Q~. However, the total charge Q = 
(Q~) + Q + is independent of this choice of phase, provided we add the appropriate quantity of the constant positive 
background charge density. 

The grand canonical potential is the relevant thermodynamic potential to describe the mechanical properties of the 
electron gas in the nanowire, and may be written 

n = -i J dED(E)ln(l + e-^ E -^). (8) 

It should be noted that Eq. (||) only holds for noninteracting electrons, since the thermal trace is taken assuming each 
fermionic mode is independent. A more general expression for Q, valid within a self-consistent mean-field treatment of 
interactions, will be presented elsewhereEJ Here it is sufficient to note that it will be argued below that the corrections 
to Eq. (&) due to charging effects are negligible. 

Eqs. (Q) and (|^) are deceptively simple, and it is worth emphasizing that /i is the asymptotic electrochemic al po itential 
of electrons injected from the reservoirs, not a local Fermi energy, as introduced by certain other authors. D(E) 
is the global energy density of eigenstates of the scattering problem, and contains all effects of multiple scattering, 
quantum interference, etc. These eigenstates are populated according to the Fermi distributions of the reservoirs. 
The occupation of a local basis of states, which are not eigenstates of H, is in general quite complicated, and will be 



discussed in detail in Sec. IV E 



In the following, it will be assumed that the volume Ld 2 of the nanowire (or the area Ld for 2D nanowires) remains 
constant under elongation, i.e., we assume an ideal plastic deformation. It can be shown that relaxation of this, 
constraint, to include e.g., a small elastic deformation, does not modify the mesoscopic effects in an essential way.Ell 
The cohesive force of the nanowire is given by the derivative of the grand canonical potential with respect to the 
elongation: 

Differentiating the above expression for £1 and performing a partial integration over E, one obtains the general result 

(10) 



F = ~ J dE f(E)lm^\ndet[S(E)]Y 



We thus have relations which give us the conductance, charge, and force as a function of the S-matrix. 



III. THE S-MATRIX OF THE WNW— GEOMETRY 



In order to obtain exact results for the tunneling and finite-size corrections to the conducting and thermodynamic 
properties of metallic nanowires, we consider the special WNW geometry illustrated in Fig. [I], for which the S-matrix 
can be calculated exactly. The system consists of a noninteracting electron gas confined by hard-wall boundary 
conditions (we will return to the question of electron-electron interactions in Sec. |VB| ). The wire cross-section 
in the 3D case is assumed to be square. A generalization to an arbitrary cross section is straightforward if the 
eigenfunctions and eigenvalues of the 2D Schrodinger emiation are known for that shape. The scattering problem for 
the 2D WNW geometry was solved by Szafer and StoneEl in connection with the problem of conductance quantization 
in 2D semiconductor quantum point contacts, and was further investigated by Weisshaar et aVEn In addition to the 
transmission coefficients calculated in Refs. ^l| and |3^, we need the reflection coefficients, i.e., the full S-matrix. The 
generalization of the method of Refs. [u 32 to 3D nanowires and to calculate the full S-matrix is described below. 

In order to calculate the elements of the S-Matrix, solutions of the Schrodinger-equation are matched at the 
transitions between wide and narrow parts of the wire. Let us regard the transition from wide to narrow first (cf. left 
part of Fig. |l|). If the z-coordinate is directed along the wire and x is an abbreviation for x in the two dimensional 
case and for (x\, X2) in the three dimensional case, describing the dimcnsion(s) perpendicular to the z-axis, the wave 
functions are given by 
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*(z <0,x) = e iK » z $ N (x) + r N ' N e- iK »' z $ N ,(x) , (11) 

N' 

*(z>0, X )=J2 t 'nNe iknZ <Pn(x) , (12) 
n 

where we have assumed an incoming wave from the left of unit amplitude, and <&jv and <fr n are transverse eigenfunctions 
in the wide and narrow parts of the wire, respectively, and are given by 

*w(aO = J-sin — (x + D/2)) N = 1,2, ... (13) 



D \ D 

/2 / 727T \ 

-smi K —(x + d/2)j n=l,2,... (14) 

in the 2D case, and by an analogous expression consisting of a product of two sine-functions for the 3D case. Here D is 
the diameter of the wide part of the wire and d the diameter of the narrow part. In the 3D case, the transverse modes 
have in general two indices, e.g., P and Q, but we can order the states according to their energy and so characterize 
them by one quantum number N; P and Q are then functions of N . The wave numbers and k n are 



K N = J- T E-— (15) 



>- 


TT 2 N 2 


D 2 


2m 

1F E - 


ir 2 n 2 
d 2 



kn = ^—E- — , (16) 

where we use the abbreviation (for the 3D case) N 2 = P{N) 2 + Q(N) 2 . 

The solution of the Schrodinger equation must obey two conditions at the transition point z = 0: 

1. Continuity of the wave function for x € [-D/2,D/2] in 2D (xi,x 2 G [-D/2,D/2] in 3D): 

& N (x) + J2 r N'N^N'(x) =^t' nN <j> n (x)G(d/2 - \x\) , (17) 
N' n 

where Q(x) is the Heavyside step function and is an abbreviation for the product 9(d/2 — \xi\)®(d/2 — \x2\) in 
the 3D case. 

2. Continuity of the first derivative of the wave function for x £ [—d/2,d/2] in 2D (xi,X2 € [— d/2,d/2] in 3D): 

K N ® N (x) - ^2r N > N K N i^ N >(x) = t' n 

) ■ (18) 

N' n 

These equations can be transformed into matrix equations for r and t' by multiplication with $>m(%) and <p m {x), 
respectively, and integration over the appropriate s-range. Using the abbreviations 

l-d/2 

PNn = / dx$ N (x)4> n (x) , (19) 
J-d/2 

Knn> — Snn'Kn , k nn i — S nn >k n , (20) 



two equations for r and t are obtained: 



l + r = pt>, (21) 
p T K - p T Kr = kt' . (22) 



Note that the p matrix is not orthogonal. The equations (|2l|) and (22) can be solved, and we obtain 



t' = 2(k + p T K P y 1 p T K , (23) 
r = pt' - 1 . (24) 



An exactly analogous calculation for an incoming wave from the right side gives 
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r' = (k + p T Kp)- 1 (k~ p T Kp) , (25) 
t = p + pr' . (26) 

The scattering matrix is obtained by normalizing the wave amplitudes with respect to current (the unitarity of S 
reflects current conservation). With 

t nN = {k n /K N ) 1/2 t nN , r NN , = (K N /K N ,) 1/2 r NN , , etc., (27) 
the scattering matrix is given by 

SW = (?}')■ ( M > 

The S-matrix for the combined WNW-geometry may be constructed from three scattering matrices U, and 
S( 2 \ describing the scattering at the WN boundary, the free propagation within the narrow section, and the scattering 
at the NW boundary, respectively (see Fig. |l|). The free propagation in the narrow section is described by the matrix 

U = ( y n ) ; X nn > = 6 nn > exp(ik n L) . (29) 



v X / 

The NW transition is associated with a matrix 

S<*> = (3°) 

which can be calculated like or can be seen by symmetry considerations. To compute the full S-matrix, we use 
the linear equations connecting the current amplitudes propagating between the individual scattering matrices (see 
Fig. [I] for an explanation of the notation) : 

b i ) = ^ 



U[%), CM., 



a t ) ~ U \ b s, 

b X \ = 5(2) ( a s 2 
b s 2 J \ a s 2 

Eliminating the unwanted variables in this set of linear equations, and rewriting the equations in the form 

%.)->{%)■ 

relating ingoing and outgoing currents, the full S-matrix is found to be 

(&n + ^^(l _ Ui 2 S < ^U2iS2 1 2)^ 1 Ui2S[ 2 1 > U2iS2 1 i 'S'i^I 1 ~ Ui2S[ 2 { } U2iS22 ) )~ 1 Ui 2 s[ 2 2 \ 

p. 
£>2l(l - ^ r 2l5'22^12»S'il^) _1 C ;r 2lS'2^ S 2 \ + (1 - 1^21 U\2 S$ ) ~ 1 U 2 1 S 2 \ U\2 S[f J 

(33) 

The operator P is a projection operator onto the undamped modes in the wide part of the wire. Note that the 
individual matrices S 1 -- 2 ^, and U are infinite-dimensional matrices describing scattering and propagation in all 

available modes, including the evanescent modes. The full S-matrix, on the other hand, connects the incoming and 
outgoing asymptotic states, and thus has a finite dimension for a given energy E, determined by the total number of 
transverse modes v with e„ < E. The inclusion of the virtual intermediate states, which describe tunneling processes, 
is crucial to solve the Schrodinger equation accurately. However, the contribution of the evanescent modes decreases 
exponentially with increasing energy. In practice, we found numerical convergence of the S matrix if roughly 20 times 
more modes than the undamped modes in the W-part were retained. 

We remark that the elements of the S-matrix can also be found by considering e.g., the transmission matrix £12 
as the sum of the directly transmitted current amplitudes and the multiply backscattered current amplitudes. This 
results in a geometric series, which can be summed to obtain the result (J33|) - 
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IV. 2D NANOWIRE 



In this section, we investigate the properties of 2D nanowires. There are several motivations to study 2D systems: 
First, a quasi-2D nanowirc could be experimentally realized in a thin metallic film. Secondly, the characteristic 
electrical and mechanical properties of a nanowire, namely, conductance quantization and force oscillations, are 
already present in 2D systems, and it is worthwhile to investigate to what extent the universality of the mesoscopic 
force oscillations predicted in Ref. |l5| depends on dimensionality. Further, 2D electronic structure is easier to visualize, 
making it simple to study the correlations between the measured quantities and the local electronic structure. Finally, 
and perhaps most importantly, certain of the phenomena studied here are directly applicable to 2D quantum point 
contacts formed in semiconductor hctcrostructurcs.Qia While the predicted mesoscopic force oscillations of order 
Ep/Xp would be many orders of magnitude smaller in doped semiconductors due to the smaller Fermi energy and 
correspondingly longer Fermi wavelength, and would likely be hidden by the much larger cohesive forces associated 
with the covalently bonded electrons of the valence band, the charge oscillations predicted to accompany the quantized 
steps in the conductance should have a comparable size in both metallic and semiconductor quantum point contacts, 
namely, of order the fundamental charge quantum e. 



A. Force and conductance 

Once the S-matrix (^) is known, the conductance and cohesive force can be calculated from Eqs. (|) and ©. 
Fig. H shows the behavior of the conductance and cohesive force as a function of the elongation of the wire. An ideal 
plastic deformation of the narrow part is assumed, which means that its area A = Ld = LqD is held constant, Lq 
being the initial length of the narrow section. The wide sections of the wire support 5 propagating modes at the 
Fermi energy; this fixes the conductance of the wire before deformation to be G = 5Gq. In Fig. g(a), the area of the 
N-part is comparatively large (3A^), while it is small in Fig. |^(b) (0.5A|.). In both cases, conductance quantization 
can be observed; whenever a channel in the N-part closes (as the diameter is decreased), the conductance decreases 
and reaches a plateau at an integer multiple of 2e 2 //i. The conductance does not show a perfect step-like structure, 
though: When the narrow part is short, there will be tunneling through the constriction before a channel opens, 
and reflection above the threshold. This will smear out the steps, leading to a rather smooth transition between the 
plateaus [see Fig. §(b)]. When the narrow part is quite long, on the other hand, there is almost no tunneling, but 
a resonant structure near the transition points occup, [Fig. ||(a)] ■ This is due to the alternating constructive and 
destructive internal reflection within the constriction.Eil 

The cohesive force is strongly correlated with the conductance. In Fig. ||(a), the modulus of the force increases along 
the conductance plateaus, while it decreases sharply at the conductance steps. The behavior is qualitatively similar 
to the result for smooth 3D geometries presented in Ref. |l5|, and to the experimental results for 3D Au nanowires. an 
Thus we see that the essential correlations of the electrical and mechanical properties of nanowires are present even 
in 2D systems, and even for abrupt geometries. For the extremely short nanowire considered in Fig. |2j(b), similar 
force oscillations correlated with the conductance steps are present (see Fig. , but they are superimposed on a much 
larger background force. The pronounced difference in the force in Figs. p](a) and |^(b) indicates a breakdown of the, 
invariance of F under a stretching of the geometry d(z) — » d(Xz), which was found within the WKB approximation,!^ 
due to strong tunneling effects in very short wires. 

The cohesive force decays to zero as the wire is elongated past the point where the last conductance channel is cut 
off, though rather more slowly than the conductance itself. F remains noticably finite even for the largest elongations 
shown in Figs. ||(a) and |](b), although G is exponentially smallH The force in this regime arises from the variation 
of the free energy due to a deformation of the geometry in the region classically forbidden to electrons, and is clearly 
larger in the shorter wire [Fig. §(b)], where tunneling effects are more important. This effect is simple to understand: 
even when the probability to tunnel all the way through the narrow section is exponentially small, the probability to 
enter the narrow section need not be small [see Fig. ||(d)], so the electron gas is still sensitive to its shape. 

It is clear from Eq. ( JTo| ) that all states with energy smaller than the Fermi energy contribute to the total force. On 
the other hand, the graphs show that the force oscillations are correlated to the behavior of the conductance, and 
thus must be due to states near the Fermi energy. In Fig. |](a), one can see that even the resonant structure in the 
conductance is reflected in the force, leading to sudden changes in its derivative. 

Fig. |i| shows the conductance and cohesive force as a function of elongation for wires with different outer diameters. 
While the curves are distinct at the beginning of the elongation (where inner and outer parts have comparable 
diameters), there is almost no difference at higher elongations, and this is not only valid for the conductance^ but 
also for the force. This shows that even a narrow constriction coupled to contacts with infinite width can be accurately 
modeled by wires of the type we regard here. 
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B. Mesoscopic force and charge fluctuations 



In order to understand the overall behavior of the cohesive force, it is useful to consider a semiclassical expansion!!; 
of Eq. (||). Formally, the density of states may be written D{E) = dN(E)/dE, where N(E) is the total number of 
states with energy less than E in the system. The behavior of N(E) for 2D domains jwift hard-wall boundary 
conditions was first investigated by Weyl,E-3 and was further developecL-by Kac and others.E3E3 The Weyl expansion 
of N(E) for a 2D domain with a polygonal boundary with n corners istj 

A dA 

N(E) = —k 2 E -—k E + 2C + SN(E), (34) 

ZTT Ztt 



where k E = y 2mE/h 2 is the wavevector associated with energy E, A is the area of the domain, dA the circumference 
of the domain, and C is a constant depending on the topology, in this case 

n 1 2 

C = Y,^f^> (35) 

where ji is the interior angle of corner i jr and SN(E) is a fluctuating term associated with the discreteness of the level 
spectrum, whose energy average is zero£3 Note that Eq. (|34|) includes a factor of 2 for spin. The shift in the total 
number of modes due to the presence of the sharp corners is given by C = 1/9 for the WNW geometry of Fig. 0. 
Integrating Eq. (|^) by parts and taking the limit T — > 0, one finds 

Ef *t,t*s ,™ ivEp A , 2E F . 

F 



O = - / N{E) dE = — -fA + -r^-dA - 2E F C + SCI, (36) 

A ip 3A p 



where SQ is a fluctuating mesoscopic correction. Differentiating with respect to L, subject to the constraint 
A = const., yields a semiclassical expansion for the force 

where SF = —d(5Q)/dL. The leading order term in the semiclassical expansion of the force is the surface tension, 
F sur f . For the WNW geometry, the circumference of the nanowire is dA = 2L + 2(D — d) + const., and one obtains 

^frH)- < 38) 

This indicates that the surface tension increases with increasing conductance (the Sharvin formula reads G/Gq ~ 
2d/\F in 2D) and with decreasing area of the wire. The surface tension is plotted in Fig. || as a dashed curve. The 
exact force oscillates around it. 

The Weyl expansion for the electronic charge of the nanowire for T — > is 

2ttA OA 



X F A 



F 



(Q~) = -eN(E F ) = -e[-^- — + 2C)+ 5Q . (39) 



In Fig. ||, the force fluctuations SF, calculated by subtracting the surface tension from the exact force, and the charge 
fluctuations SQo, calculated by subtracting the term in parentheses in Eq. ( |39"| ) from the exact charge, determined 
from Eq. are shown for the nanowires considered in Fig. [| The scale of the force oscillations is E f /^f (see Fig. 
|j) , similar to the result for smooth 3D geometries presented in Ref. [l5|. Strongly correlated with the force oscillations 
are charge oscillations of order the fundamental charge q uantu m e. The charge fluctuations SQo are calculated in the 
absence of screening. Screening will be considered in Sec. [VB| ; here it suffices to note that the screening properties of 
the 2D electron gas in e.g. GaAs are quite poor, so that the charge fluctuations in small-conductance 2D quantum 
point contacts should be essentially unscreened. The predicted charge oscillations should be experimentally observable 
with a local probe, such as a single-electron transistor. 
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C. Topological force 



A close examination of Fig. |2j indicates that the exact force deviates significantly from the surface tension for very 
small elongations: The discrepancy is roughly 3Ef/\f in Fig. ||(a) and XAEp/Xp in Fig. ||(b) for L — > Lq. This 
behavior is to be contrasted with the results for smooth geometries, presented in Ref. Gil in which the mesoscopic de- 
viations from the semiclassical result were always found to be < Ep/Xp (hence the term universal) . The nonuniversal 
corrections to the force in the WNW geometry at small deformations have a topological origin: Before deformation, 
the perfect wire has a smooth boundary, but as the wire is stretched, eight sharp corners develop. The sharp corners 
lead to a shift in the grand canonical potential 

Ar! top = -2E F C = ~E P . (40) 

However, the electrons incident from the reservoirs can only resolve the individual corners when they are separated 
by a distance greater than or of order A^. Thus one may expect the topological correction to f2 due to the sharp 
corners to evolve smoothly with the initial deformation, and to saturate when the corners become separated by an 
amount of order A^ . In order to see this explicitly, let us consider A£! t0 p as the work done by a topological force 
-Flop — F — F sux f which is present for small deformations L « Lq. For infinitessimal deformations, the adiabatic 
approximation becomes exact, and gives F(Lq) — — 4Ep/3\p. Thus F top (Lo) = AEpD /3XfLq. Estimating the work 
done by F top to be -F top (L )AL/2, and equating this to the change in the grand canonical potential Af2 t0 p, gives 
AL = XpLo/iD as the deformation regime where the topological force is important. The corresponding change in 
diameter is 

Ad = -DAL/L = —Xf/3, (41) 

indicating that the shift in the free energy of the system due to the introduction of sharp corners indeed saturates 
when the separation between the corners becomes comparable to the Fermi wavelength. That the force associated 
with this change of topology can be large compared to Ep/Xf is a remarkable result. 



D. Comparison with the adiabatic approximation 

In previous -theoretical investigations of the conductance and cohesive force in metallic nanowires, an adiabatic 
approximation^ was employed. The WNW geometry clearly violates the conditions of validity of the adiabatic 
approximation; it is nonetheless instructive to compare our exact results to those obtained within an adiabatic 
approximation for this abrupt geomety, in order to evaluate the importance of interchannel scattering. 

The solution of the 2D Schrbdinger equation can be written formally as ^(x, z) = (f> z (x)ip(z). The adiabatic approxi- 
mation consists of neglecting the derivatives d<fi z (x)/dz and d 2 (j> z (x)/dz 2 (which would be justified if \dD(z)/dz\ <C 1), 
so that the Schrbdinger equation decouples into separate transverse and longitudinale wave equations, 

-|^C(z) = ^)C(z), (42) 

-^^(z) = (E-E u (z))tP(z), (43) 

The 2D scattering problem then decouples into independent one-dimensional scattering problems for each channel, 
and the scattering matrix reduces to a 2 x 2 matrix for each channel, which can be computed e.g., via the WKB 
approximation.!^ It is then straightforward to calculate conductance and cohesive force using the formalism described 
in Sec. || 

While the adiabatic approximation should be a good approximation for boundaries with smoothly varying diameter, 
this condition is certainly not fulfilled in the WNW-geometry. Employing it nonetheless, the equation of motion for 
the longitudinal coordinate becomes just the ID Schrbdinger equation for a square barrier: 

(-T-3TS + 9 (*M L - Z ^ V ») ^ = ^( z ) (44) 
\ 2m oz z J 

with V v = EpZj£- — jjy) and E v = E—Ef ^"d 2 ' ^ e transmission and reflection coefficients are simply calculated 
using the continuity of the wave function and its derivative at the potential steps. We find 







(l-e 2 *^) (kl-Kl) 



tu = 



-k 2 + e 2lk - L k 2 - 2k v K v - 2e 2ik » L k„K v - K 2 + e 2ik * L K 2 ' 
-4a ih » L k v K v 

-kl + e 2lk » L kl - 2k v K v - 2e 2ik - L k v K v - K 2 + e 2lk » L K 2 , ' 



(45) 



where K v — k F \J E v / E F and k v — k F J (E v - V V )/E F . 

Inserting the S-matrix elements (|4^) into Eqs. (Q) and (|Io| ), one obtains the adiabatic approximation for the 
conductance and force. These are shown as dotted curves in Fig. ||. The adiabatic approximation captures some of 
the qualitative features of the exact solution (solid curves), but is not quantitatively correct. Since the conductance 
is more or less quantized, the discrepancy with respect to the exact solution has to be small. Not so for the force; 
the adiabatic approximation clearly fails to describe correctly even the leading order contribution to the force, the 
surface tension, when several channels are transmitted, especially for very short nanowires. Interestingly, the adiabatic 
approximation overestimates the effects of tunneling: both the conductance plateaus and the force oscillations are 
better defined in the exact calculation than in the adiabatic approximation. Perhaps the most striking conclusion 
that one should draw from Fig. |^ is that even for the worst-case scenario of an abrupt geometry, the adiabatic 
approximation works remarkably well for nanowires of length L > Af. 



E. Local density of states 



The correlations between the quantized steps in the conductance and the oscillations of order E F /X F in the force 
were interpreted in Ref. |l5| in terms of a simple physicaL-oicture, which was essentially the converse of the conven- 
tional interpretation. The conventional interpretationlErEil of the experiments of Refs. ^,|| is that the jumps in the 
conductance are due to abrupt changes of the structure of the nanowire at the atomic level, e.g., through the breaking 
of bonds, and that these structural rearrangements manifested themselves as abrupt changes in the cohesive force. 
While certainly a plausible viewpoint, the strong statistical evidenceo'LTuI for conductance quantization has no natural 
explanation within this framework. In order to substantiate the converse point of view, that the conductance channels 
themselves can be interpreted as mesoscopic bonds providing the cohesion, it is worthwhile investigating the local 
electronic structure of a nanowire within the jellium model. 

The electrical conductance is determined by the electronic structure of the nanowire in the vicinity of the Fermi 
energy. While the total force and charge clearly depend on all the states with energy below E F [c./., Eqs. (^) and 
(fTo|)], the mesoscopic force and charge oscillations, because they are correlated with the conductance steps, must 
also be essentially determined by the electronic states at the Fermi energy. The spatial character of the electronic 
wavefunctions associated with scattering states of a given energy E is naturally represented through the local density 
of states (LDOS), D(E,x,z). The LDOS of 2D quantum point contacts with smooth boundaries has recently been 
investigated by Ulreich and ZwergerO 

While Refs. p"8|-p0| regarded the density of states to be independent of the z-coordinate and to depend only on the 
x-coordinate in the narrow part of the wire, this is certainly not the case for the full solution of the problem. The 
z-dependence of the LDOS will be especially important when not only the width but also the length of the nanowire 
are on the nanometer scale, as in the experimentally relevant geometry. Since we know the exact S-matrix and the 
eigenfunctions for the WNW geometry, we can calculate the LDOS. The LDOS is obtained as a sum of the densities 
created by the different incoming channels. Only reflected or transmitted waves from the same channel will superpose 
coherently. In the wide section on the left side of the constriction, the LDOS at the Fermi energy is given by 



D(E,x,z<0) = tJ2 



N=l 



oc 



+ 



oo 



-iK N ,z 



N 



N' = l 



N 



N'=l 



y/hKjqi /m 

2 



(46) 



where Sy are the submatrices of the S-matrix 



and we have chosen z = at the boundary between wide and 
narrow parts. For the calculation of the LDOS in the narrow part of the wire, we need the current amplitudes bg and 
bg 2 in the narrow section as a function of the incoming current amplitudes Og and a^ 2 (see Fig. [j]) . In the system of 
linear equations (|3l|), the unwanted variables have to be eliminated, and the remaining equations may be rewritten 
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We find 



J s 2 



t = 



'Si _ 



= t 



b s 2 



til ti2 
til £22 



with the components of t given by the known matrices (p8|), S^ 2 " 1 (30) and U ( |29| ) as 

-1 



(2) 

12 ' 



tn = (l - S ( 2 ^Ui2S[ 2 1 > U2ij S^i , 

tl2 = S^U^ ^1 — U2lS^2 Ul2^j S 

t 2 l = ( 1 — S[iU 2 lS2 2 Ui2 ) S[iU2lS2i , 



^22 — ( 1 — S[i U21S22 U12) S12 



(47) 
(48) 

(49) 
(50) 
(51) 
(52) 



The LDOS in the narrow part of the wire is obtained as 



JV ma: 



D(E,x,0 < z < L) = - 



N=l 



E 



(tll)n' 



N 



\/hk n </i 



--<t>n>(x) + (t 2 l) n , 



N 



j — ik n r {z — L) 



■<f>n'(x) 



E 



(*12) r 



/ z 



'N 



e ~ik n , (z—L) 
= <j> n >{x) + (t 22 ) n , N — u , , 4>n'{x) 



\/Tik n , Jr 



(53) 



The LDOS is of course symmetric about the axis z = L/2. 

Fig. U shows the LDOS at the Fermi energy for 4 different elongations of a wire with a relatively long N-part (it has 
the same parameters as the wire in Fig. ^](a), so one can compare the conductance and force at these elongations with 
the LDOS). It can be seen immediately that the LDOS exhibits a highly nontrivial structure. The number of maxima 
in x-direction in the N-part of the wire reflects the number of open channels transmitted through the constriction; 
in Fig. ||(a) there are three open channels (and thus G » 3Go see fig. ||), in (b) there are two open channels and in 
(c) only one channel is left. Fig. ||(d) shows the exponential damping of the wave function just after the last channel 
has closed. So the highest open mode dominates the transverse structure in the LDOS in the narrow part of the 
wire; this can be understood considering that normalization of the wavefunctions to unit current implies that the 
wavefunctions are proportional to 1/ \fd — d n , where d n is the wire diameter at which the nth channel opens. Note 
that even in the tunneling regime [Fig. |B|(d)], the probability for an electron to enter the classically forbidden region 
can be non-negligible. 

As the conductance is a property of states at the Fermi energy, we should expect not only the number of transmitted 
channels to be reflected in the LDOS, but also the resonant structure exhibited by the conductance in Fig. ||(a). This 
is indeed the case; when we are at conductance maxima, the LDOS inside the N-part is much larger than that at the 
minima, and is very strongly modulated in the z-direction. This is because constructive interference of the multiply 
backscattered waves leads to a quasi-bound standing-wave state at the conductance maxima, while the conductance 
minima are associated with a condition of destructive interference [compare Fig. |^(b) (conductance minimum) and Fig. 
^(c) (conductance maximum)]. For the first conductance maximum after the channel opens as one widens the narrow 
section, there is one maximum in the LDOS of the N-part in the longitudinal direction; for the nth conductance 
maximum of the resonant structure, there are n longitudinal maxima of the LDOS. We thus see 5 maxima in Fig. 

. n . . . 

The electronic structure of the nanowire shown in Fig. r\ gives vivid pictorial support to the claim advanced in Ref. 
|l5| that conductance channels should be interpreted as mesoscopic bonds, which provide the cohesion of the system. 
The claim advanced here that the electronic structure in such a nanowire is dominated by quantum-confinement 
effects rather than by atomistic effects is in agreement with STM studies of electron "corrals" on Cu surfaces SB 



V. 3D NANOWIRE 



While the results for 2D nanowires presented in the preceding section are interesting both in illustrating the 
generality of the mesoscopic phenomena in question and for their relevance to experiments QHrpoint contacts in quasi- 
2D electron gases, the only experiments to date on the mechanical properties of nanowiresiB involve 3D metals. In 
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this section, we consider a 3D wire with WNW-geometry and square cross-section. For a square cross section, many 
modes are doubly degenerate, as in the case of cylindrical symmetry, leading to conductance steps of both 2e 2 /h 
and Ae 2 /h. It would be possible to lift this degeneracy by considering a wire with a rectangular rather than square 
cross-section. The formalism to compute the S-matrix (cf. Sec. Ill) and to obtain the force, charge, and conductance 
(see Sec. ||) is the same as in the 2D case, although one needs to include more evanescent modes for an accurate 
computation of the S-matrix than in the 2D case. It is straightforward to extend the present calculation to wires 
of arbitrary cross-section, if the cigenfunctions and eigenvalues of the 2D Schrodinger equation are known for that 
cross-section. 

Although the exactly solvable geometry considered here is somewhat special due to the presence of sharp edges, 
the gross behavior of the conductance and force is, similar to that observed experimetally in 3D metallic nanowirestrtj 
and calculated for smooth, adiabatic geometries The agreement of the present results for the WNW geometry with 
the experimental results of Ref. || is poorer than for the smooth geometries considered previously,ll3 indicating that 
the experimental geometry is undoubtedly much smoother than that considered here. 



A. Force and Conductance 

Fig. U shows the conductance and tensile force as a function of elongation for two model 3D wires. Under elongation, 
the narrow section is assumed to deform plastically, i.e., its volume V = Ld 2 = L D 2 is held constant, where L is 
the initial length of the narrow section. The inclusion of an additional, small elastic deformation can be shown not 
to modify the mesoscopic effects in an essential way.E3 The comparison of the exact results shown here to the results 
of an adiabatic approximation is similar to that in the 2D case (cf. Fig. ^), so for clarity we have not shown them for 
the 3D case. 

In Fig. ^|(a), a nanowire with a volume of 4A'p is shown, while a shorter wire with a volume of 1.25A|n is shown in 
Fig. [|(b). The width of the wide sections is D = 1.76A^ , which fixes the number of asymptotic propagating modes 
to be 6. In Fig. ^(a), one sees conductance plateaus at G = 1,3,4, 6 x Go, with a pronounced resonant structure 
superimposed due to multiple reflection at the abrupt junctions between wide and narrow sections. The sequence 
of degeneracies corresponds to the square symmetry of the cross-section (cylindrical symmetry, on the other hand, 
gives G = 1,3, 5,6, ••• x Go, see Ref. [l5]). Just as in the 2D case discussed above, the force exhibits mesoscopic 
oscillations of order Ep/Xp, which are correlated with the conductance steps: \F\ increases along the conductance 
plateaus, and drops abruptly at the conductance steps. The resonant structure in the conductance is also reflected 
in the first derivative of the force, particularly on the last conductance plateau. The force is similar in magnitude to 
that calculated for a smooth geometry in Ref. 15 and observed experimentally in Au nanowires in Ref. ^| the forces 



required to cut off the last two conductance eigenmodes are of order 1.25Ef/ Xf and 2.5Ef/Xf, respectively (recall 
that Ep/Xp « 1.7nN in Au). 

The surface tension has been plotted fpc .comparison as a dashed curve in Fig. ^. It has been computed analogously 
to the 2D case from the Weyl expansiorJi§Ej of the grand canonical potential, 

1 a _ o 

{l - E '\^ F v -W F dv + Ix-F^ CiLi]+ ' m 

p r edges i 

where V is the volume and dV the surface area of the nanowire; the topological terms are proportional to the lengths 
of the edges Li, and the appropriate constants are Ci = 1/4 for an edge with an inner angle of 7r/2 and Gj = —5/36 
for an edge with an inner angle of 37r/2. The surface tension, or semiclassical approximation to the force, is obtained 



from the derivative of the semiclassical approximation to [the term in parentheses in Eq. (54)] with respect to L, 
which yields 

E F ( nd nd 2 2 . 
F surf = _ + __---__ (55) 



Xf \2Xf 2LXf 3 2.7L 
for a 3D nanowire with WNW geomet ry an d square cross-section. Aside from the initial deformation, where the 



topological force is important (cf. Sec. IV B ), the force exhibits oscillations centered about the semiclassical result 
(dashed curve). 

In Fig. |^(b), a shorter wire whose conductance versus elongation matches the experimental curve shown in Fig. 1 
of Ref. is shown; Lq was chosen such that the elongation required to decrease the conductance from 6Go to is 
2Xf ~ lnm. We see that the conductance and force are correlated in a similar way. Due to the shorter length of the 
constriction, the conductance steps are smeared out by tunneling and by above-threshold reflection; the plateau at 
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G = 4Go and the associated structure in the force are no longer visible. The resonant structure in the conductance 
is also suppressed, except on the last plateau, where the narrow section is longest. The overall magnitude of the 
force is larger than for longer constrictions, due to the increased surface tension, and the total elongation required to 
break the nanowire is less. (jNpte that the effective surface tension can be reduced by up to a factor of 5 by including 
a small elastic deformation. However, the oscillations of the force around the semiclassical approximation are of 
the same order as in the previous case (see Fig. 0). The pronounced difference in Figs. ||(a) and ||(b) indicates a 
breakdown of the invariance of F under a stretching of the geometry d{z) — > d(Xz), which was found within the WKB 
approximation,E_3 due to strong tunneling effects in very short wires. 

As in the 2D case, the force decays to zero with increasing elon gation after the last conductance channel is cut 



off, although more slowly than does the conductance itself (see Sec. [V A for a discussion). F , 
even for the largest elongations shown in Fig. ^, when the conductance is exponentially small.l 



emains non-negligible 
3 A similar effect was 



observed experimentally (c.f. Fig. 1 of Ref. m, although it is not clear whether the effect was above the noise level. 



B. Charge oscillations and screening 

The charge on such a 3D nanowire may be calculated from Eq. (0), as in the 2D case. The charge on the nanowire 
changes as the wire is elongated due to surface terms and mesoscopic oscillations. The Weyl expansion for the 
electronic charge of the nanowire is 




The term in parentheses in Eq. (56) varies smoothly as the geometry of the wire is altered, while SQq describes the 
mesoscopic oscillations associated with the opening or closing of discrete transverse modes. 

In Fig. 0, the mesoscopic charge oscillations SQq, calculated by subtracting the term in parentheses in Eq. (pq ) 
from the exact charge computed via Eq. (ffl), and the force oscillations, calculated by subtracting the surface tension 
from the total force, are plotted for both wires shown in Fig. |as a function of elongation. As in the 2D case, there is 
a strong correlation between the two quantities, and the charge oscillations are of order the fundamental quantum of 
charge e. The force oscillations are, as in the case of an adiabatic geometry studied in Ref. [l5|, of order Ep/Xp, aside 
from the nonuniversal topologic al cor rection occuring for small deformations from an ideal wire, which was discussed 



in detail for the 2D case in Sec. IV C 



In an interacting system, screening of the charge oscillations will occur. The net charge SQ including screening can 
be estimated within the Thomas-Fermi approximation as follows: 

SQ = SQ - e 2 D(E F )SV, (57) 

where SQq are the charge oscillations in the noninteracting case regarded above, D(E) is the density of states (in- 
tegrated over the length of the constriction) and the potential SV due to the charge imbalance on the wire can be 
estimated as 

SV=% (58) 

We have introduced a phenomenological quantity C corresponding to the total capacitance of the inner part of the 
wire to its surroundings. Eqs. fl57j ) and ( |58"| ) can be used to compute self-consistently the charge on the wire within 
linear response: 

sq ^ + Jd°e fV c - < 59 > 

The charge in the noninteracting case SQq has already been computed and discussed (see Fig. Q). 

On dimensional grounds, the capacitance obeys C = aL, where a is a geometrical constant of order 1 (a may 
depend logarithmically on the ratio D/L). The density of states can be computed by a spatial integral of the LDOS 
over the narrow section, or from the asymptotic scattering phase shifts via Eq. (JsJ) - The later definition includes the 
contribution of the Friedel oscillations induced in the wide sections. In Fig. ||, both densities of states are shown 
for wires with the same parameters as above. They are approximately equal, indicating that the excess charge SQ 
induced on the nanowire under deformation resides mainly on the narrow section. The intricate resonant structure 



13 



in D(Ep) occurs due to the formation of quasi-bound states due to multiple reflection at the junctions of the wide 
and narrow sections (see Fig. ||), and would not be present for a smooth geometry, such as that studied in Ref. [lj]. 
Aside from this resonant structure, the overall magnitude of D(Ep) can also be determined from a Weyl expansion. 
To leading order, one finds D(Ep) ~ (AL / EpXp)G / Gq. Inserting this expression into Eq. (|59|), one finds 

(60) 



^ 1 + 0.66r s G/G a' 

were r s is the dimensionless electron gas parameter, which takes values between 2 and 6 in metals. This indicates that 
the screening of the mesoscopic charge fluctuations is poor for wires with small conductance. The screened mesoscopic 
charge fluctuations should be measurable with a local probe, such as a single-electron transistor. 

The screening of the predicted charge fluctuations should be even weaker in 2D GaAs quantum point contacts, due 
to the large dielectric constant, which enhances the capacitance. In the 2D case, the constant 0.66 is replaced by 0.52 
in the denominator of Eq. (|60|), and the geometrical factor a — > ea, where e s» 13. This indicates that the predicted 
charge fluctuations are essentially unscreened for small-conductance QPCs in GaAs. 

Let us finally add a comment on the effect of screened electron-electron interactions on the free energy. Within 
linear response, the Coulomb energy associated with the mesoscopic charge imbalance is given by 

An c = - J d 3 x D ■ E. (61) 

In the discrete potential model introduced above to treat screening, this leads to 

^ = crSwr (62) 



The details of the derivation of Eq. (|62j) will be given elsewhere.ti] Equation (|62|) indicates that the mesoscopic charge 
fluctuations of order e lead to a negligible correction to the free energy of the system, even in the limit C —* of 
perfect screening, justifying the independent-electron model of nanocohesion. 



VI. CONCLUSIONS 



In the present article, we have investigated the conducting and thermodynamic (including mechanical) properties 
of metallic nanowires with a wide-narrow-wide geometry, using a free-electron model with hard-wall boundary con- 
ditions. All properties of the nanowire were related to the electronic scattering matrix, which was evaluated exactly, 
including all effects of tunneling and interchannel scattering. The present results confirm the central conclusion of 
Rcf. which was based on an evaluation of the scattering matrix within the adiabatic and WKB approximations, 
that closing a conductance channel by stretching a metallic nanowire requires a force of order Ep/Xp, or roughly a 
nanonewton in monovalent metals, independent of the total number of conducting channels. 

In contrast to this "universal" behavior under a smooth deformation of the geometry, we have shown that the 
force associated with a change in topology can be large compared to Ep/Xp, and indeed comparable to the total 
(macroscopic) cohesive force. 

In addition, we predict that the net charge on a nanowire exhibits oscillations of order the fundamental charge 
quantum e, which are synchronized with the force oscillations and conductance steps. These charge oscillations 
should also be present in quasi-2D quantum point contacts in GaAs heterostructures, and should be experimentally 
detectable using a local probe, such as a single-electron transistor. 

A final word should be added by way of addressing the central controversy of this fieldp- wiiich. can be stated 
as follows: Is it the atomic structure (i.e., the bonds) which determines the conduct ancejEnEfE] or should the 
conductance channels themselves be thought of as mesoscopic bonds which provide the cohesion, and thus determine 
the structure? ■— ■ 

The conductance channels are the eigenstates of the electronic scattering problemjlj and are thus the appropriate 
states to describe both the transport and thermodynamic properties of a nanowire, which is an open quantum 
mechanical system. These scattering states are linear combinations of local bonding states, so there is no fundamental 
contradiction between the two viewpoints stated above: one is free to look at the problem in a localized basis of bonds 
or in a basis of extended electronic eigenstates. 

However, since an exact solution of the many-body Schrodinger equation for a nanowire is beyond our current 
capabilities, one is forced to. make certain approximations, which are convenient in the basis of choice. Thus molec- 
ular dynamics simulationsEilcJ typically neglect any quantum coherence between different bonds, and amount to a 
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computational version of the classical ball-and-stick model of atomic structure, where bonds are described by a short- 
ranged empirical inter-atomic potential. This is an uncontrolled approximation, which should be adequate to describe 
covalent bonds in an insulator, but its applicability for monovalent metals with nearly spherical Fermi surfaces like 
Au and Na is questionable. 

The molecular dynamics simulations involve empirically determined short-range interatomic potentials whose char- 
acteristic length and energy scales mimic the quantum mechanics of bonding. When playing classical mechanics with 
these quantum forces, it is not too surprising if one obtains forces of the right order of magnitude. Such models 
are of course inadequate to describe electrical conduction, so to explain the observed correlations in the conductance 
and force of metallic nanawires, a quantum-mechanical model whose geometry is fit to the results of the classical 
simulation is constructedHatij The cost of working in a localized basis is thus the necessity of using different physical 
laws to describe conductance and cohesion. 

On the other hand, we have seen in the present article (see also Ref. |l5|) that the observed correlations in the 
conducting and mechanical properties of metallic nanowires can be accounted for naturally in a single quantum 
mechanical model, which treats the mechanical and electrical properties of the system on an equal footing. In order 
to solve the quantum scattering problem, we have neglected the discrete atomic structure, working in a jcllium-like 
model, which is equivalent to assuming that the only effect of the lattice is to modify the electron's effective mass. 
This should be a rather good approximation for simple metals like Na and adequate for noble metals like Au. Defects 
in the atomic structure of the wire or roughness, in its surface introduce additional scattering, which can also be 
included in the jellium model in a natural way.Ej A drawback of the jellium model, or at least of the assumption 
employed here and in Ref. [l^ that the positive background deforms continuously as the nanowire is elongated, is the 
inability to describe the hysteretic behavior found in the experiment of Ref. 0. The claim advanced here that the 
electronic structure in such a nanowire is dominated by quantum-confinement effects rather than by atomistic effects 
is in agreement with STM studies of electron "corrals" on Cu surfaces.cS 

In the end, the merits of the jellium model vis a vis an atomistic description must be decided based on its predictive 
power. An interesting prediction of the jellium model discussed above is the existence of mesoscopic charge fluctuations 
of order e, which are strongly correlated with the force fluctuations, and synchronized with the conductance steps. 
These charge fluctuations are a collective effect, which can not be described in a ball-and-stick picture of bonding. 
It is incumbent on the experimenter to verify or falsify this clear prediction of the jellium model. 
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FIG. 1. Schematic diagram of the WNW-Geometry. In the upper part, the wire geometry is sketched. D and d are the 
diameters of the wide and narrow parts of the wire, respectively. Under elongation, it is assumed that the area of the narrow 
part of the wire is conserved, Ld — LoD = A, where Lo is the initial length. In the 3D case, the wire is assumed to have a 
square cross-section, and the volume of the narrow part is held constant during elongation, Ld 2 = LoD 2 — V. The lower part 
of the figure shows the scattering scheme: scattering matrices (S^ 1 ' and S 1 ' 2 ' at WN-transitions and U for the narrow part 
of the wire) characterize the transmission and reflection of current amplitudes denoted as arrows in the diagram. The total 

and a~c 



S-matrix relates the outgoing current amplitudes 6j and b $2 



to the incoming current amplitudes a,uu u, S2 . 
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FIG. 2. Electrical conductance and tensile force for two different 2D wires with WNW-geometry as a function of the length 
L of the narrow part. Assuming area conservation, L is varied from a perfect wire (where narrow and wide parts have the same 
diameter) until the last conductance channel breaks, (a) A wire with a relatively large area of the narrow part (A = 3.0\%) 
and (b) a wire with a smaller area (A = 0.5Af.) are shown. The diameter of the wide part is 2.9Af, fixing the total number 
of asymptotically propagating modes to be 5. The dotted curves show the adiabatic approximation to conductance and force, 
the dashed curves give the surface tension, the leading order contribution in a semiclassical expansion of the force. The arrows 
indicate the geometries used in Fig. H 
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FIG. 3. Electrical conductance and cohesive force for 2D wires with various outer diameters D. The area of the narrow part 
was held constant, so the initial length (when D = d) of the wires is different. The area of the narrow part is 3.0Af- as in Fig. 
^|(a). The curves are vertically offset. 
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FIG. 4. Force and charge oscillations for two different 2D WNW wires as a function of elongation. Force oscillations are 
shown in the upper and charge oscillations in the lower half. The wire parameters are the same as in the preceding figures, (a) 
shows a wire with a larger and (b) a wire with a smaller area of the narrow part. 



20 





x [Ap] x [A F ] 



FIG. 5. Local density of states D(Ef, x, z) at the Fermi energy in a 2D wire. The area of the N-part and the diameter of 
the W-part have been chosen as in Fig. |^(a) (A — 3.0A F and D — 2.9Af) to make comparison with Fig. ^|(a) possible. The 
diameter of the N-part is (a) 2.4Af, (b) 1.9Af, (c) 1.4Af, (d) 0.9Af- The corresponding elongations are marked as arrows in 
the upper part of Fig. ^(a). Black areas correspond to D = 0.0/(Ef\%), white areas to D > 1.0/(EfX f ), contours are drawn 
at equally spaced values of D between these two limits. 
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FIG. 6. Electrical conductance and tensile force for two different 3D wires with WNW-geometry and square cross-section 
as a function of the length L of the narrow part. The width of the narrow part is determined by a constant -volume constraint 
Ld 2 = LqD 2 — V — const, where Lo is the initial length of the constriction before deformation. In (a) the volume of the 
narrow part is given by V — 4Af- while it is smaller (V = 1.25Afr) in (b). The dimensions in (b) are comparable with those of 
the Au wire studied experimentally in Ref. 4. The dashed curves give the force expected by the surface term plus topological 
correction. 
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